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SUMMARY 


This paper presents a general theory for small deformations 
of anisotropic elastic shells subjected to arbitrary temperature 
distribution. 

The shells are assumed to be homogeneous through the thick- 
ness except for the coefficient of thermal expansion, but are 
otherwise unrestricted in homogeneity and isotropy. 

The theory is specialized for arbitrary cylinders, circular 
cones, spheres, and axially symmetrical orthotropic shells of 
revolution. In addition, two special problems are analyzed. 


INTRODUCTION 


The behavior of structural elements subjected to heating has recently 
assumed importance in many fields. Although most materials now in use 
are at least approximately isotropic in nature, there are instances where 
more general analyses based on considerations of anisotropy and non- 
homogeneity are necessary. Materials often become anisotropic in manu- 
facturing processes such as cold rolling or stretching. In addition, 
temperature gradients may cause the elastic constants to vary, thus 
introducing non-homogeneity. 

The theory presented in this paper is restricted to small deflections 
of thin elastic shells. The Kirchhoff assumption, normals to the middle 
surface remain straight, normal, and inextensional during the deforma- 
tion, is used throughout. Although they are assumed to be constant 
through the thickness, the elastic coefficients of the material are arbi- 
trary functions of the surface coordinates x and y of the shell. On the 
other hand, the thermal coefficients and the temperature distribution are 
unrestricted. 


NOTATIONS 


X, Y,Z = Rectangular coordinates. 
t,3,& = Unit vectors along the No aed, 
axes. 
7 = Position vector of the point X, 
Naas 
x,y,z = Shell coordinates based on the 
lines of principal curvature. See 
Section IB. The symbols x and 
y are also used as arbitrary curvi- 
linear coordinates in Section I A. 
s = Arc length. 
EH, F,G = Coefficients in the ‘first funda- 
mental form’’ of a surface. See 
Eqs. (3) and (4). 
S = Surface area. 
nm = Unit normal vector to a surface 
at any point. 
-e,f,9 = Coefficients in the “second funda- 
mental form” of a surface. See 
Eqs. (8) and (9). 
"1,72 = Principal radii of curvature de- 
fined by Eq. (11). 
M = Mean curvature of a surface de- 
fined by Eq. (10). 
K = Gaussian curvature of a surface 
defined by Eq. (10). 
A, B = Positive functions of x and y de- 
fined by Eq. (13). 
a, B, y = Lamé coefficients for orthogonal 
curvilinear coordinates. See Eqs. 
(19) and (20). The symbol a is 
also used as the apex angle of a 
cone in Section III B, but no con- 
fusion should result. 
h = Thickness of the shell. 


Tr, Ty, 4 


= Stress components at point (2, y, 
Tay) Tyzy Tzx 


z). See Fig. 1. 
Bye cy LV ye 
Q:, Q,; = Tensions, shears, twisting mo- 
i Ne wes IM « ments, and bending moments de- 


fined by Eq. (25). 


P,, Py, P, = Components of resultant external 
force per unit area on middle sur- 
face of the shell. 


kz, Ry = Components of the external cou- 
ple per unit area on the middle 

surface of the shell. 
Ex, €y, é,| 


= Strain components at point (z, 
Yay) Vyz) Vex 


2) 


U,V, W = Components of the displacement 
vector in the directions of the co- 
ordinate lines. 


u,v, w = Components of the displacement 
vector of the middle surface in 
the directions of the coordinate 
lines. 

Pi, Qi, Ri 
P2, Qo, Rey = Coefficients defined by Eq. (85). 
Sts du ME 
T = Temperature measured above an 
arbitrary zero. 
b;; = Elastic coefficients. See Eq. (86). 
c; = Thermal coefficients. See Eq. 
(36). 
Ci, Cz = Constants defined by Eq. (38). 
U, = Strain energy density. 
= Strain energy of the shell. 
Function defined by Eq. (44). 
= Shell coordinate used in place of y 
for the cone and the shell of revo- 
lution. See Sections III B and D. 
s,c,t = Abbreviations for sin a, cos a, and 
tan a respectively in the equa- 
tions of Section III B. The sym- 
bols s and c are also used as 
abbreviations for sin x and ctn x 
in Section III C. 
r = Coordinate used for shell of revo- 
lution. See Fig. 5. 
Vi, V2 = Alternative set of displacement 
components for shell of revolu- 
tion. See Fig. 6. 
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1. GEOMETRICAL PRELIMINARIES 


The theory of shells presented here is based 
upon the differential geometry of surfaces and 
elasticity theory. The geometrical topics discussed 
in this chapter are presented in books such as those 
by Struik™* or Graustein.® The following treat- 
ment is essentially the same as that presented by 


Langhaar and Boresi.® 


A. GEOMETRY OF A SURFACE 


A point in space may be located with reference 
to a rectangular cartesian coordinate system (X, 
Y, Z). The radius vector 7 from the origin to the 
point may be written as follows: 


FHiX+jV+kZ (1) 


where 7, j, & are unit vectors directed along the 
positive X, Y, Z axes respectively. 

A surface in space is represented by the equa- 
hones xX —X(z, y), Y=Y(2, y), Z=Z(a, y), where 
2, y are arbitrary curvilinear coordinates. Then a 
point on the surface is denoted by the relation 
=F (texts) 

The square of the infinitesimal distance ds be- 
tween two neighboring points of the surface is 


gs — ar dr = (f,dr +f, dy)? (2) 


(3) 


where subscripts x and y denote partial derivatives. 
By comparison of Eqs. (2) and (3), we obtain 


or 
ds? = EF dx? + 2F dx dy + G dy’. 


I 
31 


2° ls = Oe as Y/+ 2,7 
Bre, 7.2, 
yi ly = Xx, a MG ra 2, 


(4) 


on a> ans) 
I 
3 


ll 
= 


Equation (3) is called the ‘‘first fundamental form”’ 
of the surface.” Equation (4) shows that # and 
G are positive. Since ds? is positive, Eqs. (3) and 
(4) yield HG—F?>0. 

An infinitesimal element of area on the surface 
is given by the formula 


* Superscripts in parentheses refer to references listed 
in the Bibliography. 


dS =/ EG — F? dz dy. (5) 


Hence, the area of any part of the surface is given 
by the relation 


S = | [VEG — FP dedy. 


The unit vector 7 normal to the surface is given 
by the formula 


(6) 


r Treg SS They ia SSR 


at chile AOE 


(7) 


This equation fixes the positive sense of 7% as well 
as its direction. 

Another important relation in the theory of 
surfaces is the following: 


—drF-dh = edx? + 2fdxdy + g dy’. (8) 


Equation (8) is known as the “‘second fundamental 
form” of a surface. The coefficients e, f, and g are 
given by the following formulas: 


1 Spee | 
aa 
VAT TNS a 
: ne RZ 
I= iG S Z a @ 
ii Xy Yuy Ly 


The extreme values of the curvatures of normal 
sections of the surface at a given point, called the 
‘‘principal curvatures” of the surface, are denoted 
by 1/ri and 1/rz. The mean curvature M and the 
Gaussian curvature K are defined as follows: 


ie s its ce Twi? HO Wien mL) 


In the following, it is assumed that lines of principal 
curvature coincide with coordinate lines x, y. With 
this restriction, f=/=0. Then the principal cur- 
vatures are determined by the following formulas: 
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1/r, = g/G. (11) 


Equation (11) determines the signs of 7: and r2 as 
well as their magnitudes. 

A relation of differential geometry which finds 
application in shell theory is Rodrigues’ theorem.“ 
For f=F =0, it may be expressed in the form 


1/n = e/E, 


on 2 ey ie on tt Soe (12) 
Ox r OX Oy 


For orthogonal coordinates, it is convenient to 
introduce the notations 


AbaH, BY=tG. (13) 


It is shown in differential geometry that the 
functions E, G, e, and g satisfy three differential 
equations of compatibility, known as the Gauss- 
Codazzi equations. For orthogonal surface coor- 
dinates (fF =0), the Gauss equation is 


a Ay) 
eve a (14) 


_ 9 Bs) 


where K is the Gaussian curvature. The Gaussian 
curvature is a bending invariant. For f=F=0 
(that is, for coordinate lines which coincide with 
lines of principal curvature), the Codazzi equa- 


tions are 
apa ) Ps eee a ( B ) _ 1 oB 
Oy ( 1 ony YT? Oy } Ox 1 am val Ox . (15) 


B. GEOMETRY OF A SHELL 


It was noted in Section A that 7=7(2, y) repre- 
sents a surface in space. In addition to coordinates 
x, y, a third coordinate z is required to represent 
a shell. Thus, a point in a shell may be represented 
by the relation 7=7(a, y, z). If coordinates 2, y, 2 
are orthogonal at every point, they are called 
orthogonal curvilinear coordinates. Then, since 
the vector derivatives 7,, 7,, 7, are tangent to their 
respective coordinate lines, 

7, Fy =F) +f, = FF, = 0. 


(16) 


The square of the infinitesimal distance be- 
tween two neighboring points is 


ds? = dr. dF = ¢,.dz +7, dy +f,d2)> “QD 
With Eq. (16), Eq. (17) may be written simply 
ds? = 7, °F, dx? +7, > fy dy? + Ta: fede a) (18) 


Let the following notations be introduced: 


Om me Tah he, Biz tpety, vy? = f, fe 1(19) 
With this notation, Eq. (18) becomes 
ds? = 0? dz? + 6? dy? + 7? dz’. (20) 


The coefficients a, 8, y are functions of a, y, z. 
They are called ‘‘Lamé coefficients.” 

A special type of orthogonal curvilinear coor- 
dinates is employed in shell theory. Let x and y 
denote surface coordinates on the middle surface 
of the shell. Let z denote the distance of a point 
measured from the middle surface in the direction 
of the normal to the middle surface. Positive z is 
measured in the positive sense of the unit normal 
n of the middle surface (see Eq. 7). The lateral 
boundary of the shell is z= +h/2, where h is the 
thickness of the shell. In general, h is a function 
of x and y. When z, y, and z are defined in this 
way, they are called shell coordinates. 

If the shell coordinates are orthogonal, the co- 
ordinate lines on the middle surface coincide with 
the lines of principal curvature. This follows from 
a theorem of Dupin which states: The surfaces of a 
triply orthogonal system intersect on their lines of 
principal curvature. For orthogonal shell coordi- 
nates, the square of the infinitesimal distance ds 
between two neighboring points of the middle 


surface is 


Gs = Age 2b ay (21) 


For orthogonal shell coordinates, the Lamé co- 
efficients (Eq. 19) become 


a A(l f= 2/11), Pi Bil a z/T2), om L (22) 


where 7; and rz are the principal radii of curvature 
of the middle surface (see Eq. 11). 
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A. EQUILIBRIUM RELATIONS 


The stress notations used are shown in Fig. 1. 
The stress c, 1s normal to a plane which is perpen- 
dicular to the x-axis, and the stresses T2,, Tz are 
tangent to this plane and directed in the y and z 
. directions, respectively. It is shown in the theory 
of elasticity™ that the shearing stresses obey the 
- following relations: 


Try Tys, Tae = Taz, Tae Tay: (23) 


Love™ has derived the differential equations of 
equilibrium for any orthogonal coordinates. For 
shell coordinates (y=1) in the absence of body 
forces Love’s equations are 


0 fe) 0 
a (Box) + ey, (aray) + Ee (aT 22) 


Oa 0a op ee 
i a Pea ee 
te) 0 (6) 
ax (BT xy) saa (ay) ee (aBTy2) 
x Oy Oz 
fey 
0p ESD Pega OC wg 
Sige UT Oo aa: Te Fats g 
fe) 0 te) 
ae (Br22) + by. (atys) + ae (aBo;) 
0a op x 
sy 52 
Zz 
Oz 
Py 
y 
we 
+ 
Figure 1 
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These equations are exact only if 2, y, z are orthog- 
onal shell coordinates for the deformed shell. For 
small displacements, the effects of deformation on 
the equilibrium relations are usually neglected. 

Thetensions Nz, N,, shears. Nz,,..Noz, Qe, Qs, 
twisting moments M,,, M,.:, and bending moments 
M,., M, per unit length of the middle surface 
may be expressed in terms of the stress com- 
ponents o,, Gy, Tzy, Tzz) Tye by considering an in- 
finitesimal element of the shell (Fig. 2). The 
positive senses of Nz, Ny, Nay, Nyz, Qz, Qy are 
indicated in Fig. 2, where double-headed arrows 
denote moments with positive sense given by the 
right-hand-rule convention. 

In Fig. 2, the total tensile force on the differ- 
ential element in the x direction is N,B dy, where 
N, is the tension per unit length on the plane 
perpendicular to the z-axis. Also, with the nota- 
tion of Fig. 1, the total tensile force may be written 
in the form 


[os 8 dy de = dy [ Bo. dz. 


Hence, 
h/2 
N= = [6 a, dz = fle (1 — 2/72) dz. 


Similarly, Ny, Ney,\..-, Q, may be expressed. 
Thus, we obtain the following equations: 


Ne = [oo (1 — 2/re) dz, (25) 


Figure 2 
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h/2 
Ny =| oy (1 — 2/r1) dz, 
—h/2 


h/2 
ioe 


—h/2 


Na — 2/12) dz, 


h/2 
Nye =[ Tay (1 = 2/11) dz, 
—h/2 


h/2 
o,2 (1 — 2/re) dz, 


—h/2 


M, = 

(25) 
h/2 

M, ide oye (1 — z/r) dz, 


h/2 
My =[ try 2(1 — 2/12) daz, 


h/2 


n/2 
M yx =[ Tr 2 (1 — 2/ri) dz, 
h/2 


Q; Bile (1 — 2/re) dz, 


h/2 
h/2 

ae =| Tyz (1 — 2/11) dz. 
—h/2 


These representations of tensions, shears, and mo- 
ments are similar to those used by Fliigge. 

The equilibrium equations for a shell are de- 
rived by Langhaar and Boresi.“ The results are 
repeated here. 


aA aB 
-*_ (BN.) + y SANs) + Na ay Nee 
= PO. JAR tet) (26) 
ry 
P 9 aB aA 
“5p AON ea) ee oy AN gor Nas go ae nae 
aes AARP ast (27) 
Sony, LARP f= 0 (28) 
aA 
ane 7 (BM:) age (A Myz) =: Mig =e oy 
aB 
a meyer (29) 
au 
a 7 (BM) geese 7 (AMA) ae 
_M, 05 aa ABR,=0 (30) 
eee alae 
hoe Ob g Hele ae (31) 


In these equations P,, P,, P, denote the compo- 
nents of the resultant external force per unit area 
on the middle surface, and R, and R, denote 
components of the external couple on the middle 
surface of the shell. If R, and R, are neglected, 
the moment equilibrium equation is 


seta 
dx | A 


[-2- BM.) + F- (AM) 
+ Ma — Mapp | 
+ ay Ue Lee 


se 0A 
a aco mt, 


(BM) + = (AM,) 


(32) 


B. STRAIN-DISPLACEMENT RELATIONS 


Let U, V, W denote the components of the dis- 
placement vector in the directions of the coordinate 
lines. The general expressions for the strains e,, 
€y, €z) You; Yue) Yer, in terms of U, V, W, have been 
derived by Novozhilov® and Shaw.” However, 
the objective of the theory of shells is to reduce 
the problem of three-dimensional elasticity to two 
dimensions. To attain this goal, we employ the 
Kirchhoff assumption, which states that under a 
deformation, normals to the middle surface remain 
straight, normal, and inextensional. This implies 
that €,=Yy2=Yzr=0. It follows that the displace- 
ment vector U, V, W is a linear function of the 
thickness coordinate z. Also, for small deflections, 
quadratic terms in U, V, W are usually neglected. 

With these approximations, the general formu- 
las for the strain components e;, €y, Yey at any 
point in the shell in terms of the displacement 
components u, v, w of the middle surface have been 
derived by Love and by Langhaar and Boresi. 
The results are 


pei oe (vay w 
Cag + (1 — 2/r)™ eeate 
= (cael ee 
eae Ox (—) 
od See Wes yWy 
Bea (4 \+a5 Ere. +4, 


(33) 
Vy 


B, w 
oa ae ee Pan 
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—2(1 — 2z/r:)- (72 Oy aria 


ea Gps 1 r nel 

(33) 
Oe alae ae “ey, 
+ (1 -2/ny Co) 


+ ; 7 a ne 
mw 


We ve) BW, 
2 = A Bae 


—2(1 —2/r)7 ae 


For the purposes of this paper, it is convenient to 
rewrite Eq. (33) in the following form: 


bia Q:1 i zh, 
eet o/r, 1 —2/r’ 
:3 Q» D zR, 
east 1 — 2/r 1 —2/re’ 
(34) 
i a K 
at — z/r 1 — 2/r 
a 2L w zM 
1 —2/n 1 — 2/re’ 
where 
vA, Ble — ub; oo 
coe AB ry? Q: = aR Tr? 
wu 3 il i @) We 
ae A ac GE A an ec ) 
vA, AyWy 
BG) AB? 
B oa v 
eee A on UR ) 
me oo fl oe ) 
Re B oy (-)+ : : 
@ i (35) 
U zx gWy 
PaaGr, 0 AtB’ 


Brera, a): 


‘2 = eG ) (35) 


Wy 
u 4 )- 4 
A 2 So u 
a Br; Oy G ) 
BW, 
= i a a) - A B2° 


C. STRESS-STRAIN-TEMPERATURE RELATIONS 


For linearly elastic anisotropic material, the 
plane stress-strain-temperature relations for a 
shell are 


i bi; Cj Ci ih (36) 


where b;;=b; are elastic constants, and c; are 
thermal coefficients.” A repeated index is to be 
summed from one to three. In Eq. (86), the fol- 


lowing notations have been used: 
01 = Oz, 02 = Oy, 03 = Tzy) 


El 1E-5 €o = €y, €3 = Yzy- 


To obtain tensions, shears, and moments in 
terms of u, v, w, we substitute Eq. (34) into Eq. 
(36) and then substitute the result into Eq. (25). 
Thus, we obtain 


Neu) {bu [Pi +Q (1 — C,) + 1CiR] 
+ bie [P2 + Qe] + bis [K + J (1 — Ci) 


h/2 
+ Lr,C,]} =f at (1 = 2/T2) dz, 


Ny =h {boo [P> => Q» (1 a C2) a i 12C2 Ro] 
+ by [Pi + Qi] + Ues [J + K (1 =U) 


h/2 
+ Mr.C;]} — ib or tae 


Nay =h {bis [Pi sie Q: «@ a= Ci) an rCiRi] 
+ bes [P2 + Qo] + bss [K + J (1 — Ci) 


h/2 
+ LrC,]} — ‘i co (1 = 2/12) de, 


Nyc =h {bos [ot Oo (bh C2) + r2C2Ro] 
+ bis [Pi + Qi] + b33 [J + K (1 — C5) 


h/2 
+ Mr.C3}} — [ a (1 — z/r:) dz, 


a {bs [P ‘ ss 


(37) 


M, = 2 Q,Cyrire 


M, — 


My, = 


where 


Cy = (1 = r1/r2) [1 _ = In ; a 55 y 
eS 2 
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a Riri (1 i aes ) 


+ die [Pe ai Roe] + bis |=> JCyrire 


$n +n (1 


h/2 
if Tz (1 — 2/re) dz, 
—h/2 


—h3 
DOs 


+ n(s—295) 


Wy 
+ bis ae + Riri] + bog Es KCorire 


- a) |} 


h/2 
Hh oT z (1 — 2/ri) dz, 
—h/2 


— he 
STS {bs E a 


1 2rir2 C; 
a) 


+ bs. [Po + Rere] + bss [=> 


{be [Ps + + ay 2 Q.Cyrire 


+ Ir, + Mr, (1 


2 
ee Q:Cinire 
+ Riri (1 — 


h2 J Cini 9 


1 Q2riteCy 
a 


h/2 
=i c3T' z (1 — 2/re) dz, 
—h/2 
—h3 
{ba | 
~~ WT C5 6) 
h2 


oy ele DSN a chs F 


a Mro+ Lr, (1 = 


Corre 


+ Rots (1 


he Gs T1"%9 


a W2rreCs yl} 
h2 


h/2 
=I eT 2 (1 — 2/r) dz, 
—h/2 


+ Lr; + Mn(1 


(37) 


D. STRAIN ENERGY OF THE SHELL 


To obtain solutions of shell problems by energy 
methods, the strain energy of the shell in terms of 
the displacement components of the middle sur- 
face isrequired. The sum of the strain energy inte- 
gral and the potential energy of the external loads 
equals the total potential energy of the system. By 
minimization of the total potential energy, expres- 
sions for u, v, w are obtained. When uw, v, w are 
known, strain components, stress components, and 
tractions may be computed by Eqs. (34), (36), and 
(37). 
components u, v, and w are often difficult to obtain, 
approximate methods, such as the Rayleigh-Ritz 
procedure, are frequently employed. 

The expression for the strain energy of the shell 
is obtained by integrating the strain energy density 
throughout the volume. The stress components 
are related to the strain energy density” by the 
relations 


0U, 
OE; * 


oy = (39) 
Substituting Eq. (36) into Eq. (89) and inte- 
grating, we obtain 


ioe 


5 cel + GUL. 


bixeie; = (40) 
Since the arbitrary function C(7’) is immaterial in 
the application of energy methods, it is disre- 


garded. Thus we obtain 


(Gf = > bi 6€; — Gel’. (41) | 
The total strain energy is 
U = [[ [Uo a8 dx dy de. (42) 


Integration with respect to z may be performed 
readily with the assumption that b, are inde- 
pendent of z. Thus, substituting Eq. (384) into 
Eq. (41), substituting the results into Eq. (42), 


Since exact solutions for the displacement 


and integrating through the thickness, we obtain — 


a 


(Pr Qi)? +5 (P: + Qe)? 


2 f bu 
U =| fr 2 
+ Dis (Pi + Qi) (P2 +. Qe) 


+8 (K+ J)? + ba (Pi +) T+ K) 


+ bus (Ps + Q:) (J + K)} AB dr dy 
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2 
+f{s ane E ail) ie (K—rM) (Q.—rR) || AB de dy 
T1792 J 
if uaa ine ]+ Doo 2 [ (Patras) ip AB dx dy (43) 
T1719 
where 
_ bss [ (re2M +r,L)? h/2 
(Q2— Takeo)? i 2 | rite [= ayes Se Ta = 2/72) (Oil hP 2) 
3 4 er aes mL)? — ae Cis ee mM)? ] + (1 =a 2/11) (c3K =e CoQ) 
+ (1 — 2/re) (esd + €1Qs) 
<li rly) e eee) —2(1 —2/n) (@R2 + 35M) 
1/2 
—z(1 — 2/re) (e3L + c.R1)] T dz. (44) 
me [ (r,L + eM) (Pi + 1R:) see 
y rire Equation (43) is the general expression for the 
120 strain energy of an anisotropic, non-homogeneous, 
1 


Oi riBp (Q; — nk) | 


ey [ (r2M + riL) (P2 + 12h) 


T1P2 


elastic shell with arbitrary temperature distribu- 
tion. It is restricted to small deflections. For an 
isotropic shell, Eq. (43) reduces to the result 
obtained by Langhaar and Boresi.“ 


Ill. SPECIALIZATIONS OF THE GENERAL THEORY 


In this chapter, the equations developed in 
Chapter II are specialized for several common 
types of shells. 


A. CYLINDRICAL SHELL 


In this section, equations are presented for an 
arbitrary cylindrical shell oriented with its genera- 
tors parallel to the Z-axis. Then X=X(y), Y= 
Y(y), Z=2, where x is distance measured along a 
generator and y is are length around the girth of 
the cylinder. 

By Kgs. (4), (9), (11), and (13), 


Ei= tt” fF =O. 4G =a Agee 

45 
Be1, 1/rn.=0, re =r(y). Sit 

With Eq. (45), Eqs. (35) and (88) become 

12%, = Uz, Ps = Vy; Q: a 0, Qs = —w/r, 


Ry = Wee, Ry = Wy — 07/17, J = vz, 


K=yu, L=v,/r+0,, M= ws, (46) 
1+ h/ Or ) 
(Os = = _ a GPs St See some = 
1 0, (ls i i In ts h/2r Ge 


Substitution of Eq. (46) into Eq. (37) yields formu- 
las for tensions, shears, and moments: 


Ne = hibu | ue + 292 | + bref EF Fel 
== bis [ = Vz ie a (vz/7 + wy) || 
h/2 
=H aT (1 — 2/r) dz, 
—h/2 
Ny = h {bo [vy —w/r +C (w/t +rWyy —vry/7)] 
== bie Ux = bos [Vz Se Uy =i C(rWzy ax: Uy)]} 
(47) 


h/2 
= CT dz, 


—h/2 


ee =h {bss [Uy + Up = h? (vz = T Wey) /12r?| 
+ bis [Ue + A? Wee/12r] + bes [vy — w/r]} 
h/2 
Hi c3T (1 — 2/r) dz, 
—h/2 
Nae = i) {bss [v2 4 Uy + G (TWry = Uy) | oa bis Uz 
+ beg [vy — w/r-+C (w/r + ry — vry/7))} 


= [est dz, 


h/2 


hs 
Ws = 7 129 {bu [te + Wael + Die [vy + rWyy 


— vry/t| + 2bis [vs + TWr]} 
h/2 

=i «Tz (1 — z/r) dz, 
—h/2 


3 
M, = (bn (12C/h?) (ury — r?Wyy — wv) 


12 
+ bisWee + bes l oe 


h/2 
|} =! eT 2 dz, 
—h/2 


(Tuy — 1? Way) 


M oy ae a {2bs3 [ve Hrwoy] = bis [Uc tren] 
+P bes [Dy + PWyy = Ore /T]} 
h]2 
=i c3T'2z (1 — 2/r) dz, 
—h/2 
hs 
VE ae Bie: {b33 [v2/7 + Way + (12Cr/h?) (uy 


— TWey)| + bis Wee + be3 (12C/h?) (vy 


h/2 
— Pg — ae Ug) =a c3T'z dz. 
—h/2 


20 


(47) 
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By Eqs. (26) to (31) and Eq. (45), the equilibrium 
equations are 


ON.  ONvs a 
a =P ay + P,=0 
ON wy ON, a Qy ee 
Ox ap oy r Bs Caan 
& aif dQ, + Ny + yee = 0 
ot Oy r 
aM. aM. ey 
£2 Ue = 
pot a, ee +R = 0 
oe Hi 
ot oy SE get ae 
tN ae ry IN sg = te . 


iso, Eqs. (32) and (45) yield 


eM, OM y- OM wy 
Ox? Oxdoy Oxdoy 
7M, Ny, i 
hye aNeeeetane feat GO) 


The strain components, by Eqs. (84) and (46), are 
€g = Uz — 2Wee, 


er 
Yy —w/ (12) ~~ (Wy — oF y/7), 


ey 
(50) 
Yoy = V2 + TUy/(r ee OL TP Wey) 


z( af Jw 
r—z we 


By Eqs. (43) and (46), the strain energy for the 
cylindrical shell is 


U =[[zh {bir Us? + Doe (vy — w/r)? + bas (02 


+ Uy)? + 2bive (Vy — w/rT) + Qists (Vz 
++ Uy) + 2bes (vy — w/r) (uy + v2)} dx dy 


3 
a. {bu (We a 2U2Wz2/7) 


— Boy (12C/h?) (w/r + rwy — vr,/r)? 
+ bas [3 (2/7 + Way)? — (12C fh?) (uy 
— Way)? + bie [Weady/? + Wee (Way 
— vry/1?)| + 2b13 [(uz/1) (¥2/7 + Wey) 
| + Wee (Ve/7 + Wry)] + Qbes [(v2/7 


+ Way) (Vy/7 + Wyy — vr, /r?) 
+ 120 /1*) (ty = rWy) (Wy + w/P 


— vr, /r)]} dx dy Hf fd dx dy, (51) 


where 


h/2 
L= ral [((1 — 2/r) (cit, + Cody + C302) 
—h/2 


+ catty — cow/r — 2(1 — 2/r) (csWay 
ats C302/7 =| C1W2z) me (CoWyy a C3Way 


— cvr,/r?)| T dz. (52) 


B. CONICAL SHELL OF CIRCULAR CROSS SECTION 


The middle surface of the shell is given by the 
equations X =z cos a, Y=z sin a cos 6, Z=x sin 
a sin 6, where x, 6 are shell coordinates indicated 
in Fig. 3. 


Figure 3 


By Eas. (4), (9), (11), and (13), 


HE=1, F=0, G=(as)?, A 


1/r; = 0, 


a | 
(53) 


where s=sin a, c=cos a, t=tan a. With Kq. (53), 
Eqs. (35) and (38) become 


Drees re = at =1, 


5 
Ve 
Pi = tz, P2 =—, 
xs 
U WwW 
= 0 Po eee ae 
Qi ae x wee 
Wee Wz 
Ry aia Wray Ry 7TY x28? + a , 
v Ue 
J=%v,—-—, K= \ 
ae ne (54) 
v v Wao We 
L a i 2 ae = 20? 
aoe et xs es 
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ee bala C, = 0, 
xs xs (54) 
$ at 1 + h/2at ) a 
Cae dee he 5 ae 


Substituting Eq. (54) into Eq. (37), we obtain 
tensions, shears, and moments as follows: 


h2 u 

Vie h bu [us a sar Y= na ES 
+4 4] + bis [ve — +t 
xs at By as 


gfe ae Nh v Wx6 5 Wa \|} 
c 12270? (v. x a5 Cc xe / JJ 


h/2 
=i! aT (1 — 2/at) dz, 
—h/2 


Ne=h [bia Ue + doe |= = (1 
ws 


-o(2--$)+e(m +) 


Vv Ue 
+ bys [ve -2 +4 (1 — 0) 


( _ we }}- pe 
eee ib (cat Ue 


oe [ nes ] Ee 
Nig =h | bis Wee 122t Wea + bog xs 


h/2 
sth c3T' (1 — 2/at) dz, 
—h/2 


Ve h {bs ee Ee 4h 
zs 


=) 2) £10 (we + = 


x 4b 


)] 


Ue 


Vv 
+ bis [ve —2 + (1 = 0) 


—( meee c 
a wet - | fot 


h3 
M,= say {hn ae ar LEW an! 


+ bu |e + ¢( un, 28) 


(55) 


+ 2b [ve — + ms (wa = a (55) 


h/2 
rh aT 2 (1 — 2/2t) dz, 
—h/2 


hs 12éC w 
Bas {be he [uF 
4 = tw. | == bis Wer 
+ bys | (0. — 2 + Be — Se 
i x G LCG) 
12Ct (ue _—s TW We )}} 
ts h? Ss ‘t a ok 
h/2 
=i CT z dz, 
—h/2 
ae Ss | 
M.0 = 12zt jbis [Us vie rtWeel 
+ bas | + t( we + 2#-)] 
xs xs 
+ 2b5;[ 02-2 + 1(*# — 22) 
x 8 xs 
h/2 
=i csT' 2 (1 — 2/xt) dz, 
—h/2 
he 
Mo: = 5 {bi War + bog E a == 
Wee 127C. 
= ae t i Ie 


12Ct (we «FW )) 
of h? 8 Cc ef c 


h/2 
— / c3T' z dz. 
—h/2 


With Eq. (53), Eqs. (26) through (31) become 


Meg Mey MMe Sena 
oN 4 Nea fe 

races cent 
sy Be ee ee 
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OM, M, il OM ox Ms i 
ax a x a xs 06 co Q. = 0 (56) 
OM Mv. : il OM, Moz = 
0x . x os xs 06 r x Qs = 0 
Noz — Nis = ies 
at 


In Eq. (56), the terms R, and R, have been dis- 
carded. By Eqs. (82) and (53), the moment 
equilibrium equation is 


7M, 1s 2 OM, 4 | te (Moz + Me) 
Ox? ie ake as 0x00 
= 1 OM, | il 0 (M 26 a= Moz) 
La Or as 06 
ie 07, Ne he 
F 2252 ae? a= t 2 fae = 0. (57) 
By Eqs. (84) and (54), the strains are 
Ey = Uz — 2 Wee | 
eee u—w/t - 2 ( Wee ) 
me ae a xr—z/t x—z/t Wer oe 
beanie: 0 ue/S z ( (58) 
OD aan x—z2/t ie 
vs Wee We ) us (100 We )) 
c XC x—z2/t £ 


sy, 


The strain energy for the conical shell is obtained 
by substitution of Eq. (54) into Eqs. (48) and 


(44) as 
1 2 
emote) 


5 = ={[= [bu ust + 


+ ba (0e— 2 + HL) + 2h 


+u- 2) + ris ue( v2 — ates 


(ot 
thee eee! +f Ge 24 i |b & 
sk, Suse | + bn WBC 


— E Pani! 
at h? ag at 


2 & 
a t(w. ae ) + bss pe [ (w- 


XS 


+ 2bes : (u 
x 


Sie ts)? | + 2bie 


Wea 
| we +e + 28] 
40 085 xst 
+ 2bi3 
a6 


aG us) (,,-2 
t Wee + xt es a 


ees UR a 
s 
= tw.) (ug — 2tW29 + tw») |} xs dx dé 


+ [ [12s dx dé (59) 
where 


h/2 2 Va 
[= Pf 16 1- =) (cme + 2% 
h/2 at xs 
Ue u 0) 
+ Cz — C3 ~) + eg + ¢ ae C2 a 02 
Zz Wee Zz We 
aa C2 Wi A C3 Ua eae 
L ws” “Ss 4 0} 
v 
—2(1 — z/at) | ciWer + 63( 02 — = 


hee we.) : |} re. 


é 106 at 


(60) 


C. SPHERICAL SHELL 
The middle surface of a spherical shell with 
center at the origin is given by the equations 
X= Hsmnzcosy, Y= Rsin x sin y,| 
(61) 
Z = Reoosx | 
where FR is the radius of the middle surface, 2 is 


the colatitude, and y is the longitude (see Fig. 4). 
By Eqs. (4), (9), (11), (18), and (61), 


hah 1G. = bh sinta, f = 0, 
(62) 
A=R, B=Rsinz, n=n= — 
With Eq. (62), Eqs. (35) and (38) become 
a, gle ay EM en ae 
Py a R ) P» Rs F) Q: R p) 
uc +w Wee (63) 
OQ» a R ) Ri = R? ) 
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W. WC = OC Uae Oe h/2 
R, = => + alr oie (63) = esT' (1 + 2/R) dz, 
—h/2 


Res? Re 
_ Uy Com Oe Way — WyC h3 
“ets Rs’ Bee Nie 2 ies ! M, = 2 soar (hu (Uz — Wee) + bie ( 
(12R? 
Way — CWy — Uy 
MOS or Rigia || ees a + cu — cw. — MH) + bra (ve (64) 


pape QZwWy + Uy — 2Way »} 


s 
h/2 

| Tz (1 + 2/R) dz, 
—h/2 


hs 


= 12k? soar {be Ca Wen) + bee (t+ 


G 


26Wy + Uy — 2Wey )f 


s 


+ cu — cwWz — 


—ct+ 


h/2 
Sil coTz (1 + 2/R) dz, 
—h/2 


Meiditt= os a8 
zy = uz = TOR? | if is = eg) 


nS 
Figure 4 es bag (4+ + cu — cwz — Woy ) 
s & 
where s=sin z, c= ctn x. Substituting Eq. (63) + bss (v.-co+ Zetia ty ey ») 
into Eq. (37), we obtain tensions, shears, and 
; Be h/2 
moments as follows: —| caTz (1 + 2/R) dz. 
—h/2 
ei ( v 
N; == 3 -- mats 5 
R you (ust wv) + Op 8 ary Substitution of Eq. (62) into Eqs. (26) to (31) 
af yields the following equilibrium equations: 
si w) + bis ( te cv)| 
s ON; in ON 
Na ox = s ; ay GUVs Ny) 
~ fot A+ e/R) de, 40, + RP; = 0, 
h ON, id 
Ny = [biz (te + 00) + bre (2% + on oer I Ny + 2¢N xy 
+w) +da( +». - ev)| hig ides, 
OQ; 1 0 
h/2 80s 4 1 Oe _ i, + ™,) 
at @T (1 + 2/R) dz, | Y (65) 
a + ¢cQ, + RP, = 0, 
a ae 
Nay = Nox = (bus (Ue + a) + bos oe te +o (M. — M,) 
L 8 Oy 


64 
+cu+w + bss a Ey en \ ( ) — kQ. + KR, = 0, 
S J 
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(65) 
— RQ, — RR, = 0. 


By Eqs. (82) and (62), the moment equilibrium 
equation is 


eM, L mM, 2 May 
Ox? ie oy? + s dxdy ap? race 
= M,) +2272 + (M, — M.) 
—R(N,+ N,) + @P, = 0. (66) 
Equations (34) and (63) yield 
ee ws 2Wee 
omer | Rte R(R-+.2)’ 
meric Vy 
me OR 4 Rs a R = Z 
ce, — ve Uy 
elk Rs 
22 


(pa oy (haz) (Way — CWy). 


The strain energy for the spherical shell is ob- 
tained from Eqs. (43), (44), and (63) as 


v= ={{= {ou (w Se Hee + bee 


2 
2,5) 


+w+cu) 


+w+eu) Sf ak 


Vy 


Bi,’ wv) ( : 


+ 2bis (Uz + w) ( a0 he cv) 


+ 2s (2% + w + cu) (“4 +0, 


— cv)| sdx dy + [sie {bu (uz Wee)? 


+ boo (- cee teanc ihe ee “ty 


2 


Behe ( Uy + 2ewy — Way iy ov) 


Ss 


+ 2bi2 (Ue — Wee) (4 +eu-cw. ae ) 


Une We Oy, == Deg 


s 


— 2b13 (GH = Wire) ( 


+o, —- cv) + 2bo3 (» + cu — cw, 


Wy )( Uy + 2cwy — 2Wey 
s 
Se \ ele oy 2 
BCE )| sdx dy +f f{IR S dz dy, 


where 
a ea eee | al 
[= ae 1 R CiUr + : 


;, + 2/P) 


R [ eve +c (w+cu) + ¢3 (. 
Uy z2(1+2/R) 
gt s ) fi? 


Ee 
w 
+ 6 
8 


) 
ae Ce cu) + ¢3 (ve — 0; 


(68) 


2cWy 


: =— )] T dz. 


For isotropic shells the results derived in Sections 
A, B, and C reduce to those obtained by Langhaar 
and Boresi.@ ® 


QWWey — 


+ (69) 


D. AXIALLY SYMMETRICAL ORTHOTROPIC 
SHELL OF REVOLUTION 


The shell of revolution is shown in Fig. 5. The 
middle surface of the shell is described by r=r(x) 


ve 


ao 


Figure 5 
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and Z=Z(x). A point in the surface is located by is brora( Urs 4 Tate )} | 
the coordinates x and 6 where « is arc length from rry r (74) 


the origin measured along a generator to the point 
and @ is the angle between the X-Z plane and the 
plane containing the point and the Z-axis. From 
Fig. 5, 

X =rcos8, 


Z=Z. (70) 


Equations (4), (9), (11), (13), and (70) yield 


Y=resin@, 


E=1, F=0)G="7, 4A =—1, 8 =97; 


(71) 
L/t, = %iLee = Tink, Lite Za t 


Consider an orthotropic shell of revolution 
which undergoes axially symmetrical deformation. 
Then 


fa) 
v= ay = bis = bas = C3 = Noy = Nys 


= M, = M,, = Q, = 0. (72) 
With these restrictions, Eq. (35) becomes 
. ) 
Pi =z, P. =0, ha 
ry 
ae ue aa y Ras — +) + ee 
; "2 (73) 
Rs a ~( a= Ss ratte) 
, ry 
J=K=L=M =0. J 
Substitution of Eq. (73) into Eq. (37) yields 
N; = h {bu (uw. = = (1 — Cy 
+ mG Ee aa (—) i wee |) 
eye 
7 tip) 
h/2 
aa eT (1 — 2/12) dz, 
—h/2 
Kee h {ba la c1) ( =: —- = ) 
=i nCo(+ = + tw) + bie [ ue 
Hint Vig 
Ww h/2 
= = || —f et (el me 2/11) dz 
hs 120 
M, md —5— {Ou ([u. at h2 : wr | (74) 


tn [2 (1) toe] [1 22222) 


h/2 
a Tz (1 — 2/re) dz, 
—h/2 


hs 12 Ei 
—~s5 {bn Es Corie ( = =. 


12r, r2 


nf 4a) (1 aE 


+ bis E sie Ue (—) =f ves] 


h/2 
ee eTz (1 — 2/r) dz. 
—h/2 


With Eq. (71), the equilibrium equations (Kqs. 
(26) to (31)) become 


te) ip 
oy (Ne) — Fag — = Ge eee 


0 is r 
ws 2s ue = 75 
aq (0s) +7 Na +—- Net oP, = 07s ae 


2 GM.) = Mota 70 


Substitution of Eq. (71) into Eq. (82) yields 


pie M, 
Ox? 


at Te ie 7 (ee — Me) = ig lox Cie Me) 


f= Ne NG PP =, (76) 
Hit r2 


The strains, given by Eqs. (34) and (73), are 


€z = Uz — 


bese) teed 


1— 2/1 
7 (2 Ul x 2.) (77) 
6 = _— 
12 1- 
Lf aie 
ye + ratte) 1— Te 
Va) 0. 


By Eqs. (43) and (73), the strain energy of the 


shell is 
ape | 
T1 2D if Ue) 


offal 
= al rdvdé 


q 
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[ wetrin (-) + ree | 


T1192 
Peed, ae w 


[2 —n2(2)—nm]) 
boo li) Urs Tae Ne 
aR 5 a 
1 Hn fet T 
12C, ( Ur; W PQUry Toler y] 
= - — — = Wr 
h? r ee oe Sag) r 


| we erate x () ae ra | 


a ( 
ry 


[we + Te |) dnd + ff Irdedo (78) 
where 


< =f. (a = 2/73) (1 — 2/72) i, 
eel —| 2/71) a( a ~) 


T 
yi) (-a~) —22 (1 
— 2/7r1) (i + rive) —z(1 


— 2/2) C1 [ u (—-) + wee || T dz. 


Equations (74) through (79) may be expressed 
in terms of displacement components V; and V2 
in the axial and radial directions, respectively, as 


+ [fs ie 


(79) 


Figure 6 


selected by Reissner.“® Referring to Fig. 6, we 
obtain the relations for u and w in terms of V; and 
V2 as follows: 


U 


Visiny + V.cos ¥| 


(80 
Vicosy — Vesin y| 


WwW 


where y (x) is the angle between the normal to the 
middle surface and the Z-axis. In addition, the 
radii of curvature may be expressed in the form 
ye Tee BLD 
ri Ve ek tah 


(81) 
It is convenient to introduce the quantities 


e = (Vi)z cosy — (V2). sin y| 


82 
¢ = (Vijesiny + (V2). cos ¥| ms 


where ¢ is the rotation under deformation of the 
normal to the middle surface and « is the strain 
of a generator on the middle surface. 

With Eqs. (80), (81), and (82), the tractions 
become 


x Vo 
N,= hou Ja - — Cy) ea See | + die a) 
h/2 1 
-[ at (1 zs 2204.) a 
—h/2 r 
V C 
pi | 
=e (1 — zy) dz, 
me —h i bi # WAGE ' és) 
Ue glee. es 
(83) 
a Pweln v| + bis g cos vy} 
h/2 1 
eal az (1 = ane) dz, 
h/2 iF 
a ele WHER x ro 
a sy (be [ h? sin y Vs tan y 
su 2 | + bx | 
h/2 
4h Tz (1 — ez) dz. | 
—h/2 
Similarly, the equilibrium equations (Eqs. (75) 


and (76)) become 


a (rNz) — cosy No — r¥.Q. + 7Pz = 0, (84) 
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a (rQ,) +rvzN2 +siny No +rP.z = 0, 


2 (rM,) — cosy Mr —1Qz + 1Ro = 0, | ¢gqy 


ee: (2M, — My) 
(6) Ne 


— y,siny (M, — Ms) 
+ y.rN, +siny Ne+rP, = 0. 


The strains are obtained by substitution of 


Hqs. (80), (81), and (82) into Eq. (77). The 
results are 
Piper eee 
eles SA 
V2 — ze cosy (85) 
ee ; 
r—25iny 
Veo 0. 


The strain energy expression (Eqs. (78) and 
(79)) may be written in terms of V; and Vz in 
the form 


(Ss Ve 


Gus [ibs EEG a) bes 2 ride 


uberis Gs) 


— Ua (& bs c) | + bm [Pace ; 


= ae (= i ny y] 


+ 2bio (ee co8¥ rdx + 2m [ Irdz 
r 


where 
sot e - any.) i, 
Le lee E ( if (¢ 202) 


gi, ee ces v)| Paz. (87) 


(86) 


Solutions to problems of axially symmetrical 
shells of revolution may be obtained with either 
system of equations; i1.e., equations (74) through 
(79) or equations (83) through (87). The choice 
will depend on the type of loading and the bound- 
ary conditions that the particular problem presents. 


IV. SPECIAL PROBLEMS OF CYLINDRICAL SHELLS 


A. AXIALLY SYMMETRICAL DEFORMATION OF 
AN ORTHOTROPIC CIRCULAR CYLINDER 


Consider a circular cylinder, with radius a, 
oriented so that its axial coordinate xv and its 
circumferential coordinate y coincide with the 
principal directions of orthotropy. Let the cylin- 
der undergo axially symmetrical deformation. Then 

0 


ee aay. 


= M,, = My = 


= bi3 = bes = 63 = Nay = Nyz 


Q, = 0. 


The equilibrium equations, determined by Eqs. 
(48), (49), and (88), are 


(88) 


ON 
Ox 


dQ), 
Ox 


nM, 
Ox? 15 a 


(89) 


By Eq. (47), the tractions are 


aT (1 — z/a) dz, 


—h/2 


ll 


h | be (0-1) +bn | 


h/2 
CT dz, 
—h/2 
—heby 
12a 


(90) 
CUS SCE) 


h/2 
=f aTz (1 — z/a) dz, 
—h/2 


( 


Ey; 
12a 


—12aC bow 
h2 


a2 biya%vee ) 
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The strain energy for the axially symmetrical 
deformation of an orthotropic circular cylinder is 
obtained from Eqs. (51) and (52). It is 

| dx 


U —= [ah | Pus a bre (2 we y- 2bi. 


U2,W 


= Elta ae [bs (022 aie ietlss ) 
ae ae w? | de + 2x [alder (91) 
where 
h/2 
— m [es ( 1 = z/ ‘a) (Ux ae Zz) 
h/2 
a |r dz. (92) 


If the elastic constants b,;; are considered inde- 
pendent of x, the equilibrium equations may be 
expressed conveniently in terms of displacement 
components u and w. Substituting Eq. (90) into 
Eq. (89), we obtain 


ey ; Weer + Are — 5 Ws 
Af fee |. (aT) | (1—z/a) dz = 0, 

- nee a (C — 1) w +55 tess (93) 
M ae _ i + i vi oT de 
at + in Ee («T) | z (1—2/a) dz = 0. 


With the preceding equations, analyses of 
problems of axially symmetrical deformations of 
orthotropic circular cylinders may be obtained 
by two procedures. One may express the total 
potential energy of the shell as a function of u 
and w and then minimize this function. Alter- 
natively, one may solve the equilibrium equations 
(Eq. (93)) subject to the boundary conditions of 
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the particular problem. In the following articles, 
examples are presented which serve to illustrate 
the two procedures. 


B. CIRCULAR CYLINDER SUBJECTED TO 
TEMPERATURE WHICH VARIES 
THROUGH THE THICKNESS 


Let the cylinder be constrained at its ends so 
that there is no strain in the z direction (ie., 
u=0). This implies that the material properties, 
the radial displacement component w, and the 
temperature distribution are independent of z. 
Hence, the strain energy is 


Were 1 ht LC bsg ay ta ] 
U = anf [S ™. w? + wf eePde dx. (94) 


Since the integrand of Eq. (94) is independent of z, 
the integration is readily performed. The result is 


hbo» (1 = C) h/2 


Oe SS ale f w? + w 


| 2a 


Tae}, (95) 


—h/2 


where the length of the cylinder is 21. 

If the cylinder is subjected to a pressure P; 
on the inside and P, on the outside, the potential 
energy of the external loads is 


Q = 4ral (P; — P.) w. (96) 


The total potential energy V of the system is 
the sum of the strain energy U and the potential 
energy Q of the external loads. Hence, by Eqs. 
(95) and (96), 


V = 4rt{ Mat =O ys 
+[[) até + a(P.— P) |u| (97) 
h/2 v oO ° 


—h/Z 


By the principle of stationary potential energy, 
a necessary condition for equilibrium to exist is 


(98) 
Therefore, by Eqs. (97) and (98), we obtain 


a a h/2 


—h/2 


Substitution of Eq. (99) in Eqs. (36) and (50) 
yields the stress components as 


aby. [ h/2 
hb — CW = 2) [ etde 


+a(P:— P.)]— a, 


Oz = 


(100) 


a h/2 
= FS ee LL ete 
+ a(P; - P,) | ye 3 
0. 


Try 


Equations (99) and (100) express stress compo- 
nents oz, gy, Tzy and deflection w for arbitrary tem- 
perature distribution through the shell thickness. 


C. SEMI-INFINITE CYLINDER WITH END 
MOMENTS AND SHEARS 


Consider a semi-infinite cylinder as shown in 
Fig. 7. From Eq. (89), 


Ne _ 


ae 0. 


Qo 


“(| 


Figure 7 


Qo 


axial tension. Equation (90) then yields 
1 bie W h? 


+ : [ar (1 — 2/a) dz. 
2 


ails; (101) 


Therefore, N,=constant=O0 since there is‘no net — 


Substitution of Eq. (101) into Eq. (93) yields the 


equilibrium equation for w: 


Warner + 40°Wee + 48*w = f (x). (102) 


Here, 


2a7b116 : 


of = bre | (103) 
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The general solution of Eq. (102) is 


f (@) 
(103) 


(7) ] a 


oo [Bi cos-V/ B+ a x 
+ Brsin/ 6? + a? a] 
+ eV "= [Bs cos/ B+ o® 2 
+ Bysin~/ B+ 2 x] + F(x) J 


where Bi, B,, B3, Bs are arbitrary constants de- 
pending on the boundary conditions and F(x) is 
the particular solution corresponding to the tem- 
perature function f(2). 

Since w must remain finite, B;=B,s=0. The 
constants B, and B, are determined by the bound- 
ary conditions at +=0. 


(104) 


woo: = 0, MM, = M,, ae = (),. 


(105) 
Since B; and B, depend on F(x), consider the case 


T=T,e—*'*. Substitution of this function of 7’ into 
Eq. (103) yields 


i 
=) = re é ce 


Then, by Eq. (102), 


=) ge an C106) 
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(107) 
where 


(108) 


Equation (104) then reduces to 
mee Vs [Bi cos-~/ B+ a x 
+ By sin~/ B? + a? a] + ye-*/*. (109) 


Substituting Eq. (109) into Eq. (90), eliminating 
uz by Eq. (101), and applying the boundary condi- 
tions of Eq. (105), we obtain 


p= e GR) | 
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Substitution of Eqs. (103), (108), and (110) 
into Eq. (109) gives w as a function of x. Equa- 
tions (101) and (109) yield u,. Then Eqs. (36), 
(50), (101), and (109) give the stress components 
in the shell as functions of x and 2. 


D. NUMERICAL EXAMPLE 


In this section, numerical values of deflections 
and stresses are computed with the results of Sec- 
tion C for five combinations of elastic coefficients 
b;; and thermal coefficients c; (see Table 1), for a 
cylinder with mean diameter of 20 inches, and 
with thickness of 1 inch. For isotropic material 
(Case 1), 
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BE pes vH 


1-Yr 1 — vp?’ 


(111) 


where # is Young’s modulus, v is Poisson’s ratio, 
and k is the coefficient of thermal expansion. In 
Case 2, the thermal coefficient in the longitudinal 
direction has been reduced by a factor of one-half 
as compared to Case 1. Similarly, in Case 3, the 
thermal coefficient in the circumferential direction 
has been reduced by a factor of one-half. In Case 4, 
the ratio of circumferential stiffness to longitudinal 
stiffness has been doubled; in Case 5, it has been 
reduced by a factor of one-half compared to Case 1. 
For brevity in the following discussion, the ratio 
of circumferential stiffness to longitudinal stiffness 
will be designated simply as the “stiffness ratio.” 

Two types of loading with constant tempera- 
ture are considered: (a) a constant shear Q, is 
applied at the free end of the cylinder (Fig. 7); 
(b) a constant moment M, is applied at the free 
end of the cylinder (Fig. 7). Also, deflections and 
stresses are computed for a temperature distribu- 
tion T=T,e—*/*, where T, is a constant, x is the 
distance measured from the free end, and a is the 
mean radius of the cylinder. 

The deflection w is computed by means of Eqs. 
(103), (108), (109), and (110). Numerical results 
are presented graphically in Figs. 8a, 8b, and 8c. 
Temperature effects are not included in Cases 1, 
2, and 3 of Figs. 8a and 8b. In all cases, for end 
shear Q,, the deflection w is a maximum at the 
free end (Fig. 8a). The stiffness ratio is doubled 
in Case 4. Then the deflection w attains a value 
of approximately 57% that attained in the iso- 
tropic case. When the stiffness ratio is reduced 
by one-half (Case 5), the deflection is 81% greater 
than in the isotropic case. 

With end moment M,, the deflection w is again 
largest at the free end (Fig. 8b). However, for end 
moment, the deflection w is less sensitive to changes 
in the stiffness ratio; in Case 4, w decreases to 69% 
of its value in Case 1, and in Case 5, it is increased 
by 48% of its value in the isotropic case. 

Deflections due to temperature distribution 
T=T.e*/* are shown in Fig. 8c. As with end 
shear and end moment, the deflection w is a 
maximum at the free end. When the thermal 
coefficient in the longitudinal direction is decreased 
by a factor of one-half, compared to the isotropic 


Table 1 


Values of the Ratios of Elastic and Thermal Coefficients 
for the Five Cases of Article 14 


Case biz bee C1 S c2 A 
Novo. hy | eda” Oat eae ae 
1 a5) 1 10-5 105 
2 re: if 51076 10-5 
a a il 1075 5x1076 
4 ACs} rd 10-5 10-5 
> ca an) 10-5 1075 


case, the deflection w increases by approximately 
21% (Case 2, Fig. 8c); whereas, if the thermal 
coefficient in the circumferential direction is de- 
creased by a factor of one-half, w decreases by 
approximately 71% (Case 3, Fig. 8c). Thus, for 
the assumed temperature distribution, the deflec- 
tion w is more sensitive to changes in the thermal 
coefficient in the circumferential direction than in 
the longitudinal direction. Furthermore, a de- 
crease of the thermal coefficient in the circum- 
ferential direction results in a decrease in the 
lateral deflection w; a decrease in the thermal 
coefficient in the longitudinal direction results in 
an increase in the lateral deflection. Hence, if 
lateral deflection due to temperature is to be kept 
small, a shell should be designed so that these two 
opposing effects cancel each other. Cases 4 and 
5 of Fig. 8¢ illustrate the effect of the stiffness 
ratio. If the stiffness ratio is doubled, the deflec- 
tion w is decreased to 52% of its value in the 
isotropic case. If the stiffness ratio is reduced by 
a factor of two, w is increased to 216% of its value 
in Case 1. Hence, for the type of temperature 
distribution considered, the deflection w is very 
sensitive to the stiffness ratio. 

Numerical values of stress components oz and 
a, may be computed by means of Eqs. (36), (50), 


(101), and (109). In this example, co, and o, have 


been computed numerically for the free end (« =0) 
for all three types of loading. The numerical 
values of o, and oc, are illustrated graphically in 
Figs. 8d, 8e, and 8f. 

For end shear Q,, the longitudinal stress o, is 
zero at the free end (Fig. 8d). However, the 
circumferential stress oy, varies through the thick- 
ness as shown; it attains a maximum value on the 
inside of the cylinder. When the stiffness ratio is 
doubled (Case 4), the maximum value of oy, is 
increased by approximately 20% of its value in 
the isotropic case. If the stiffness ratio is decreased 
by a factor of two (Case 5), o, decreases approxi- 
mately 18%. Consequently, o, is somewhat less 
sensitive to variations in the stiffness ratio than 
is deflection. 


» per inch 


wb, 


Mo 


10° 


Winches per deg F 


-60 


lo 


t 
N 
9 


-1/80 


IV. SPECIAL PROBLEMS OF CYLINDRICAL SHELLS 33 


8 


(a) Radial deflection w 


due to shear vs axial 
coordinate x 


12 


x in inches 


fos 
x in inches 


(b) Radial deflection w 
due to moment vs axial 
coordinate x 
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For end moment M,, the variations of stress 
components o, and oy are illustrated in Fig. 8e. 
The axial stress component oz is the same in all 
cases. The variation of the circumferential stress 
component o, is similar to that obtained for end 
shear Q.. 

For temperature distribution T=T,e~*/", the 
circumferential stress component o, is shown in 
Fig. 8f, the axial stress component being zero 
at the free end. The largest compressive value 
of oy is attained in Case 2, its magnitude being 
approximately 21% larger than the value for 
the isotropic case. In Case 3, the maximum value 
of the compressive stress is approximately 72% 


less than for the isotropic case. If the stiffness 
ratio is doubled (Case 4), the maximum com- 
pressive stress is decreased by about 10%. How- 
ever, the maximum tensile stress increases in 
comparison to the isotropic case (Fig. 8e). If the 
stiffness ratio is decreased by a factor of two 
(Case 5), the maximum compressive stress 1s in- 
creased approximately 16%, whereas the maxi- 
mum tensile stress decreases. 

In general, the maximum thermal stress is 
small. For example, for aluminum, bn is approxi- 
mately 107 lb/in.2. Then, for the value T,=100°F, 
the maximum compressive stress at the outer sur- 
face is approximately 540 lb/in.? 


V. DISCUSSION OF RESULTS 


The geometrical preliminaries in Section I A are 
well known topics of the differential geometry of 
surfaces. Section IB explains the shell coordinates 
which are used throughout the remainder of the 
paper. The equilibrium and _ strain-displacement 
relations discussed in Sections IIA and IIB have 
been established previously in the theory of shells. 
The form in which they are presented here follows 
Langhaar and Boresi.“) 

In Section II C, a stress-strain-temperature rela- 
tion for a linearly elastic anisotropic material is 
chosen. The equations for the tractions, shears, 
bending moments, and twisting moments in terms 
of the displacement components u, v, w of the 
middle surface are developed. 

In Section II D, the strain energy of the shell in 
terms of the displacement components u, v, w is 
reduced to a surface integral by integration of the 
strain energy density through the thickness. 

In Chapter III, Sections A, B, C, and D, the 
general theory is specialized for cylindrical shells, 
conical shells of circular cross section, spherical 
shells, and axially symmetrical orthotropic shells 
of revolution. The latter portion of Section D 
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expresses the results for the shell of revolution in 
terms of axial and radial displacement components 
as selected by Reissner.°) 

As an illustration of the theory, two special 
problems for axially symmetrical orthotropic circu- 
lar cylinders are analyzed in Chapter IV. Section 
B treats an infinite cylinder with temperature 
varying through the thickness. Section C treats 
a semi-infinite cylinder subjected to end moments 
and end shears and a temperature distribution 
which varies in the axial direction. The problems 
are solved by two different procedures. The first 
problem is solved by expressing the potential 
energy of the system in terms of the displacement 
w, and then applying the principle of stationary 
potential energy. The second problem is solved 
by expressing the equilibrium relation in terms of 
displacement components uw and w. Then the re- 
sulting differential equations are solved for w. 

For the second problem (Section IVD), nu- 
merical values for stresses and deflections are 
computed. The results are illustrated graphically 
in Figs. 8a, b, c, d, e, and f for the isotropic case 
and for four types of orthotropy. 
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